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This  paper  provides  methods  and  experimental  results  for  recursively  estimating  the 
sprung  mass  of  a  vehicle  driving  on  rough  terrain.  It  presents  a  base-excitation  model  of 
vertical  ride  dynamics  which  treats  the  unsprung  vertical  accelerations,  instead  of  the 
terrain  profile,  as  the  ride  dynamics  model  input.  It  employs  recently  developed  methods 
based  on  polynomial  chaos  theory  and  on  the  maximum  likelihood  approach  to  estimate 
the  most  likely  value  of  the  vehicle  sprung  mass.  The  polynomial  chaos  estimator  is 
compared  to  benchmark  algorithms  including  recursive  least  squares,  recursive  total  least 
squares,  extended  Kalman  filtering,  and  unscented  Kalman  filtering  approaches.  The 
paper  experimentally  demonstrates  the  proposed  method.  The  results  of  the 
experimental  study  suggest  that  the  proposed  approach  provides  accurate  outputs  and  the 
proposed  method  is  less  sensitive  to  tuning  parameters  when  compared  with  the 
benchmark  algorithms. 


1  Introduction 

This  paper  provides  methods  and  experimental  results  for  recursively  estimating  the  sprung 
mass  of  a  vehicle  driving  on  rough  terrain.  An  accurate  onboard  estimate  of  vehicle  mass  is 
valuable  to  active  safety  systems  as  well  as  chassis  and  drivetrain  controllers.  These  autonomous 
systems  and  controllers  schedule  gear  shifts,  actuate  brakes,  induce  steer,  schedule  fuel  injection, 
warn  drivers  of  rollover  susceptibility,  etc.  Many  of  them  rely  on  accurate  knowledge  of  the 
mass  of  the  vehicle  to  perform  optimally.  Since  vehicle  mass  can  vary  significantly  from  one 
loading  condition  to  the  next,  the  estimate  of  vehicle  mass  needs  to  be  updated  online,  a 
constraint  which  adds  its  own  challenges  to  the  development  of  the  mass  estimator. 

A  significant  number  of  mass  estimation  algorithms  have  been  developed  for  on  road 
conditions.  The  authors  have  surveyed  this  literature  and  included  a  brief  summary  in  the 
introduction  of  a  paper  by  Kang,  Fathy,  and  Stein  (2008),  where  another  novel  algorithm  for  on 
road  vehicle  mass  estimation  was  presented.  Much  of  the  scientific  literature  for  vehicle  mass 
estimation  uses  the  longitudinal  vehicle  dynamics,  drivetrain  shuffle  dynamics,  or  combined 
lateral,  yaw,  and  roll  dynamics  to  estimate  the  vehicle  mass.  Despite  the  research  successes  for 
on  road  mass  estimation,  the  off-road  real-time  vehicle  mass  estimation  problem  remains 
relatively  unexplored,  and  existing  solutions  to  this  problem  remain  difficult  to  apply  in  practice. 
One  of  the  main  challenges  is  the  fact  that  the  motions  introduced  by  rough  terrain  are  significant 
enough  to  make  mass  estimation  based  on  longitudinal  vehicle  dynamics  infeasible;  this  rough 
terrain,  however,  makes  mass  estimation  based  on  vertical  vehicle  dynamics  becomes  much 
more  viable  due  to  the  presence  of  significant  terrain-induced  excitations.  The  overarching  goal 
of  this  paper  is  to  develop  an  accurate  and  fast  real-time  online  mass  estimator  for  vehicles 
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negotiating  rough  terrain,  and  to  demonstrate  its  viability  through  experimental  data  collected 
from  an  actual  vehicle. 

A  number  of  excellent  papers  and  patents  have  explored  the  use  of  the  ride  dynamics  for 
vehicle  mass  estimation.  Often,  these  methods  assume  that  the  terrain  profile  is  known, 
estimated,  or  measured  (Blanchard,  Sandu  and  Sandu  (2009a,b,2010),  Best  and  Gordon  (1998), 
Huh  et  al.  (2007),  Kim  and  Ro  (2000),  Lin  and  Kortum  (1992),  and  Shimp  (2008)).  Others 
assume  that  the  suspension  is  equipped  with  an  active  or  semi-active  force  actuator  which 
provides  a  known  suspension  force  (Ohsaku  and  Nakai  (2000),  Rajamani  and  Hedrick  (1995), 
and  Song,  Ahmadian,  and  Southward  (2005)).  Finally,  Tal  and  Elad  (1999)  analyzed  the 
dynamics  in  the  frequency  domain  to  reveal  important  resonance  frequencies  related  to  the  value 
of  the  sprung  mass  and  suspension  spring  constant.  Unlike  the  methods  above,  this  paper  seeks 
an  estimation  algorithm  which  does  not  require  a  priori  knowledge  or  explicit  measurement  or 
estimation  of  road  terrain  and  does  not  require  active  or  semi-active  suspension  force  actuation. 
This  study  also  aims  at  using  a  minimal/inexpensive  sensor  set  to  achieve  the  mass  estimation. 

The  proposed  solution  presented  in  this  paper  adopts  a  base  excitation  approach  to  real-time 
off-road  vehicle  mass  estimation.  It  introduces  the  base  excitation  model  which  treats  the 
vertical  accelerations  at  the  four  unsprung  masses  as  measured  inputs  to  sprung  vehicle  dynamics 
and  uses  the  governing  equations  of  these  sprung  vehicle  dynamics  for  mass  estimation.  This 
approach  furnishes  estimation  schemes  that  do  not  require  a  priori  knowledge,  measurement,  or 
estimation  of  the  terrain.  This  paper  also  explores  potential  benefits  of  using  a  polynomial  chaos 
based  approach  (Blanchard,  Sandu  and  Sandu  (2009a,b,2010))  for  recursive  parameter 
estimation  (Pence,  Fathy  and  Stein  (2010)  and  Pence,  Fathy  and  Stein  (2011))  and  compares  it 
with  benchmark  algorithms  such  as  regression  approaches  and  state-filtering  approaches. 

This  paper  is  structured  as  follows.  The  following  section  presents  the  derivation  of  the  full 
vehicle  base  excitation  model  used  in  this  work’s  proposed  estimator.  The  proposed  recursive 
polynomial  chaos,  and  benchmark  regressor  and  filtering  based  estimators  are  presented  in 
Sections  3-5,  respectively.  Section  6  presents  the  results  related  to  investigating  the  degree  to 
which  a  vehicle’s  pitch  dynamics  contribute  to  the  mass  estimation  problem  within  the  context  of 
an  estimator  that  uses  a  base-excitation  based  model.  Section  7  details  the  experimental  setup 
followed  by  comparisons  of  the  experimental  results  between  the  proposed  recursive  Polynomial 
Chaos  based  estimator  and  the  regressor  and  filtering  methods  in  Sections  8-9,  respectively. 
Concluding  remarks  are  in  Section  10. 


2  Base  Excitation  Full  Car  Ride  Model 

Vehicle  dynamics  literature  commonly  uses  a  seven  degree  of  freedom  system  to  model  the 
general  behavior  of  a  vehicle’s  ride  dynamics.  An  illustration  of  this  model  is  provided  in  Figure 
1.  The  seven  degrees  of  freedom  include  the  vertical  zs  cg ,  roll  9S,  and  pitch  cps  motion  of  the 
sprung  mass,  as  well  as  the  vertical  motions  zujb  zujr,  zu  rh  zu  rr  of  the  four  unsprung  masses. 
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The  motions  zg ji,zg fr,zg  ri,zg  rr  of  the  four  tire  surfaces  at  the  ground  are  the  inputs  to  the 
seven  degree  of  freedom  model. 

The  base-excitation  concept  treats  the  unsprung  mass  motions  instead  of  the  terrain  as  the 
model  input  (as  illustrated  in  Figure  2).  This  methodology  is  adopted  from  the  vibrations 
literature  (Inman  (2001))  and  has  been  used  for  mass  estimation  by  Pence,  Fathy  and  Stein 
(2009a  and  2009b)  and  Song,  Ahmadian  and  Southward  (2005),  but  these  methods  only 
considered  a  quarter-car  model 


LgJl 


Figure  1:  Seven  degree  of  freedom  model  of  vehicle  ride  dynamics. 


of  the  vehicle.  As  a  key  difference  from  the  traditional  full-car  ride  model,  the  full-car  base 
excitation  model  treats  the  vertical  unsprung  mass  accelerations  zujh  zujr,  zurh  zu,m  the 
longitudinal  velocity  U  in  the  x-direction  with  respect  to  body  fixed  axes,  and  the  sprung  mass 
pitch  velocity  (ps  as  measured  model  inputs.  (A  reduced  order  model  at  the  end  of  this  section 
will  require  only  the  unsprung  mass  accelerations  as  the  system  inputs).  The  “dot”  notation 
denotes  the  derivative  with  respect  to  time.  The  resulting  system  models  the  dynamics  of  the 
sprung  mass  in  the  following  three  degrees  of  freedom:  vertical  zscg,  pitch  (ps,  and  roll  6S.  The 
base  excitation  model  is  shown  in  Figure  3. 


Input=terrain  profile  Measure  with  low  cost  sensors 
Difficultto  measure 


Figure  2:  Base  excitation  concept 


Figure  3:  Base  excitation  model  of  full  car  ride  dynamics. 


The  base  excitation  model  has  a  number  of  key  benefits  for  vehicle  mass  estimation  compared 
with  the  traditional  full-car  ride  model:  the  model  avoids  requiring  knowledge  of  the  values  of 
the  unsprung  masses,  tire  stiffness  and  damping,  and,  most  importantly,  terrain  profile.  Also,  the 
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reduced  degrees  of  freedom  result  in  fewer  model  states  and  hence,  less  computational 
complexity. 

This  paper  makes  the  following  simplifying  assumptions: 

1.  Negligible  yaw  velocity  (xp  —  0) 

2.  Negligible  lateral  velocity  (V  —  0) 

3.  Small  pitch  (p  and  roll  6  angles 

4.  Linear  suspension  elements 

5.  Left-right  symmetry  in  suspension  elements,  e.g.,  the  front  left  spring  stiffness  is  the 
same  as  the  front  right  spring  stiffness,  e.g.,  kpi  —  kfr  kf. 

6.  The  e.g.  is  at  half  the  track  width  T /2,  a  known  distance  a  behind  the  front  axle  of  the 
vehicle,  and  a  distance  b  forward  from  the  rear  axle.  The  wheelbase  length  is  a  +  b. 

Applying  Euler’s  laws  of  motion  to  the  sprung  mass  of  Figure  3  results  in  the  following 
equations  which  govern  the  dynamic  behavior  of  the  base  excitation  full  car  model: 

z  —  equation :  m(zscg  —  (psU )  =  Fz 

(p  —  equation :  J^cpg  —  iVL,  (1) 

6  —  equation :  J90s  =  Me 

Here,  Fz  is  the  net  downward  force  acting  on  the  sprung  mass,  M(p  is  the  net  moment  or  torque 
about  the  y  axis,  and  Mg  is  the  moment  acting  about  the  x  axis.  These  forces  and  moments  are 
produced  via  the  suspension  elements  as  follows: 

Fz  "  X  Fi  (2) 

i=fl,fr,rl,rr 

My  •■=  —{Ffi  +  Ffr)a  +  (Fri  +  Frr)b  (3) 

( Ffl  +  Prl)  +  2  ( Ffr  +  Prr)  (4) 

The  forces  F[,  i  =  fl,  fr,  rl,  rr  are  due  to  suspension  deflections  and  velocities  and  are 
calculated  by 

T 


Pfi  •  kf  [Zu.fl  zs,cg  +  a-(ps  +  2  4"  bf 


Ffr  •  kj  I  Zujr 


Ffl  ■■=  kr 


T  \ 

zs,cg  T  acps  — '9S  J 

T 


+  bf 


7  1  zu,/r 


ju,ri 


-‘s.cg 


+  br 


b(ps  +  2ds) 

Frr  —  kr  ^ Zurr  —  Zs  cg  —  b(ps  F  F 

For  convenience,  this  section  defines  the  following  terms: 

b{zu,fl  T  zu,fr)  T  zu,rl  T  zu,rr) 


zs,cg  T  ^ 

T  •  \ 

zs,cg  T  a<Ps  —  ~2^s ) 
T  ■  \ 

zu,rl  zs,cg  —  b(ps  F  — -  9S  J 
T  .  \ 

^ u,rr  —  zs,cg  ~  b(ps  ~2^s)' 


u.cg 


<Pu  := 


2  (a  +  b) 

~{zu,fl  +  zu,fr)  +  ( zu,rl  +  z  u,rr) 

2  (a  +  b) 


(5) 

(6) 

(7) 

(8) 

(9) 

(10) 
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^  '  zs,cg  zu,cg  (11) 

®-=<Ps-  (pu  (12) 

The  terms  zucg  and  (pu  can  be  interpreted  as  the  average  vertical  position  and  average  pitch  of 

the  unsprung  masses.  Then,  the  terms  Z  and  <t>  can  be  interpreted  as  the  average  vertical  and 
angular  displacements  of  the  sprung  mass  relative  to  the  unsprung  masses. 

Given  definitions  (2)  Through  (12),  the  z  and  (p  equations  of  (1)  can  be  written  as: 
z  —  equation :  mZ  —  Fz  +  m(<pst/  —  zucg ) 
cp  -  equation:  =  Mv  -J^cpu 

The  vertical  force  Fz  and  pitch  moment  M(j)  can  also  be  written  as  functions  of  the  terms 
defined  in  (9)  through  (12): 

Fz  =  —2  (kf  +  /cr)Z  —  2  [bf  +  br)Z  +  2  (a/ey  —  h/cr)< t>  +  2  (afy  —  bbr)<t>  (14) 

Mv  —  2  (akf  —  bkr^)Z  +  2  (abf  —  bbr^)Z  —  2  (a2kf  +  h2/cr)  d> 

—  2  (a2bf  +  b2bry P 


(13) 


(15) 


Equations  (13)  through  (15)  are  independent  of  the  roll  motion  of  the  sprung  mass.  As  a 
result,  the  roll  dynamics  can  be  neglected  without  sacrificing  accuracy  in  the  calculation  of  the 
vertical  and  pitch  motion.  This  result  enables  a  lower  order  (and  less  computationally  expensive) 
estimator. 

The  state  space  representation  of  Equations  (13)  through  (15)  is  as  follows: 


-*i- 

*2 

*3 

,X4. 

(16) 


u 

—b 

u 

—b 

0 

—a 

0 

—a 

0 

zu,fl 

2  (a  +  b) 

0 

2  (a  +  b ) 

0 

2(a  +  b ) 

0 

2  (a  +  b) 

0 

1 

z u,fr 

zu,rl 

1 

1 

-1 

-1 

0 

0 

zu,rr 

.2  (a  +  b ) 

2(a  +  b ) 

2(a  +  b ) 

2  (a  +  b) 

-<PsU- 

The  states  are  x4  =  Z,  x2  =  Z,  x3  =  O,  and  x4  =  6.  The  set  of  dynamic  state  equations  in  (16) 
govern  the  motion  of  the  base  excitation  system  shown  in  (16).  These  dynamic  equations  can  be 
combined  with  the  estimation  techniques  of  this  paper  for  a  potential  solution  to  vehicle  mass 
estimation.  There  are  a  few  factors,  however,  that  motivate  reducing  this  set  of  state 
equations/sensor  inputs  in  (16)  to  a  lower  dimensional  set.  The  first  motivating  factor  is  the 
desire  to  reduce  the  computational  complexity  of  the  estimator.  This  is  important  for 
online/onboard  algorithms  which  may  have  limited  access  to  memory  and  computational 
resources  on  an  onboard  computer.  The  second  factor  is  motivated  by  the  desire  to  reduce  the 
number  of  sensors  required  to  perform  the  online  estimation.  Based  on  these  motivating  factors, 
this  section  derives  a  reduced  order  model,  but  will  also  note  that  the  lower  order  estimator  also 
has  a  lower  fidelity  than  an  estimator  which  uses  the  higher-order  set  of  state  equations  given  in 


Xi 
*2 
X 3 
X4J 


0  1 
— 2  [kf  +  kr)  — 2  (bf  + 


m 

0 


m 

0 


2  (akf  —  bk^  2  (abf  —  bb^ 


h 


]<L 


0  0 

2  (akf  —  bkr)  2  (abf  —  bby) 

m  m 

0  1 
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(16).  For  example,  simulation  results  in  Section  6  show  that  an  improved  estimate  is  obtained 
by  using  (16)  compared  to  the  reduced  order  model  when  the  vehicle  experiences  large  pitch 
values. 

To  reduce  the  set  of  state  equations/sensor  inputs  to  a  lower  dimensional  set,  this  paper  notes 
that  in  many  vehicles,  the  coupling  term  between  the  vertical  and  pitch  dynamics  may  be 
negligible.  In  Equation  (16),  the  vertical  and  pitch  dynamics  are  coupled  by  the  terms  2(akj  — 
bkr)/m ,  2(abf  —  bbr^)/m,  2 {akj  —  bk^) / ] v  and  2 {abf  —  bb^) l]v.  In  some  cases,  the 
numerator  of  each  term  is  small,  and  the  denominator  is  large.  For  example,  if  akj  —  bkr,  and 
abf  —  bbr  the  numerators  are  zero  and  there  is  no  coupling  between  the  vertical  and  pitch 
dynamics.  When  these  coupling  terms  are  not  weak,  a  stronger  assumption  is  required,  which  is 
to  assume  that  the  pitch  motion  of  the  vehicle  is  negligible.  Either  of  these  assumptions  enable 
decoupling  of  the  pitch  and  vertical  dynamics,  and  the  set  of  equations  given  in  (16)  can  be 
reduced  to  the  following  second  order  state  equations  that  govern  the  vertical  motion  of  the 
sprung  mass: 


0 

2  [kf  +  kr) 
m 


r  0  0 


L2(a  +  b)  2(a  +  b) 


0 

a 

2  (a  +  b) 


0 

a 


ZU,fl 

zu,fr 


2  (a  +  b)J 


zu,rl 

-zu,rr- 


(17) 


Here,  Bz  —  2  (bf  +  br).  The  measured  output  y  is  the  vertical  acceleration  zs  cg  of  the  sprung 
mass  at  the  center  of  gravity  location,  and  the  output  equation  is  as  follows: 


y  = 


— 2{kf  +  kr )  — Bz 


xi 

x2. 


(18) 


m  m  J 

The  known  or  measured  variables/parameters  of  equations  (17)— (1 8)  are  the  four  unsprung  mass 
accelerations,  the  distances  from  the  axles  to  the  sprung  mass  c.g.  location,  and  the  front  and  rear 
spring  stiffness  values.  The  unknowns  parameters/variables  include  the  value  of  the  vehicle 
sprung  mass  m,  the  damping  term  Bz.  and  the  state  variables  x1  and  x2. 

The  estimation  algorithms  are  discussed  in  the  following  sections. 


3  Polynomial  Chaos  and  Maximum  Likelihood 

Pence,  Fathy,  and  Stein  (2010  and  2011)  have  developed  a  recursive  maximum  likelihood 
estimator  for  state-space  systems  using  polynomial  chaos  theory.  This  estimator  builds  on  the 
polynomial  chaos  technique  for  modeling  multibody  dynamic  systems  with  uncertainties 
developed  by  Sandu,  Sandu,  and  Ahmadian  (2006a  and  2006b)  and  extends  the  batch  maximum 
a  posteriori  estimator  developed  by  Blanchard,  Sandu,  and  Sandu  (2009a,b,2010)  to  applications 
that  require  recursive  estimation.  Existing  research  has  explored  the  use  of  recursive  polynomial 
chaos  based  algorithms  for  vehicle  mass  estimation  (Pence,  Fathy,  and  Stein  (2009b)  and  Shimp 
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(2008)).  The  study  presented  here,  however,  uses  an  algorithm  that  estimates  parameters  based 
on  an  integrated  cost  function  which  is  different  from  the  previous  approaches  which  used  an 
algorithm  that  updates  parameter  estimates  based  on  an  instantaneous  cost  function.  The  earlier 
work  by  Pence,  Fathy,  and  Stein  (2009b)  showed  that  estimators  based  on  the  instantaneous  cost 
criteria  resulted  in  estimates  having  significant  variance.  Another  difference  between  this  and 
the  earlier  work  is  that  the  earlier  research  was  directed  at  quarter-car  suspension  models. 

The  polynomial  chaos  approach  is  applied  directly  to  the  system  of  Equations  (17)— (18). 
Because  the  system  of  Equations  (17)— (18)  is  linear,  the  Galerkin  approach  (discussed  in  Pence, 
Fathy,  and  Stein  (2010  and  2011))  can  be  applied.  To  apply  the  polynomial  chaos  based 
algorithms,  the  unknown  parameters  are  treated  as  random  variables.  The  polynomial  chaos 
algorithm  requires  upper  and  lower  bounds  on  the  unknown  parameters  which  can  be  inferred 
from  prior  knowledge.  This  paper  assumes  that  the  true  values  of  the  unknown  parameters  could 
potentially  be  any  values  between  the  lower  and  upper  bounds  with  equal  probability,  i.e.,  it 
assumes  that  the  prior  distributions  are  uniform.  With  known  upper  and  lower  bounds,  the 
random  variables  representing  the  unknown  mass  and  damping  terms  can  be  written  as  functions 
of  independently  identically  distributed  (iid)  polynomial  chaos  variables  %m  and  <fB,  each  of 
which  are  uniformly  distributed  over  the  interval  [-1,  1]: 

^  —  dm  +  Un^rn  HQ') 

Bz  =  dB+  vB^B  '  '  '  ' 


Here,  iim  and  i uB  are  respectively  the  known  mean  values  of  the  mass  m  and  damping  Bz 
random  variables.  The  terms  vm  and  vB  are  the  known  maximum  variation  values  of  m  and  Bz 
respectively  from  the  mean  values  \im  and  [iB ;  e.g.,  vm  =  mupperLimit  -  \im  where  mupperLimit 
is  the  largest  possible  numerical  value  of  the  mass  m.  The  goal  of  polynomial  chaos  based 
estimation  is  to  estimate  the  most  likely  realizations  of  m  and  ^B  and  hence,  by  (19),  the  most 
likely  values  of  the  unknown  mass  m  and  damping  term  Bz. 

Applying  the  Galerkin  method  of  polynomial  chaos  theory  to  the  system  of  Equations  (17)- 
(18)  (with  the  unknown  parameters  m  and  Bz  replaced  by  their  polynomial  chaos  counterparts 
given  in  (19))  results  in  the  following  set  of  deterministic  state  equations: 


X  = 


r  or 

®rxr  ^rxr 

X  + 

-1 

-A2i  A22 

-Or- 1 

u 


(20) 


Here,  0 k,k  E  {r  x  r,r,r  —  1}  is  a  matrix/vector  of  zeros  having  dimension  k,  the  identity 
matrix  Irxr  has  dimensions  r  x  r,  and  r  =  (S2  +  35  +  2)/2  where  S  is  the  user-selected 
maximum  polynomial  chaos  order,  (5  =  6  in  the  experiments  of  this  paper).  The  state  vector 
X  6  M2r  is  a  vector  of  polynomial  chaos  expansion  coefficients.  Finally,  the  matrices  A21  E 
Rrxr  and  A22  £  Mrxr  are  defined  as  follows: 


A21  A 


diag 


<(&!,()  d>2) 
(<P2,Q4>2) 


l)  2) 


(^l.Q^r) 

(®2  ,Q®r) 


,Q  = 


+  kr) 

dm  T  Vm^m 


(21) 
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A 


22  - 


Adiag 


(<*>!, GQJ  <<DlfG<D2>  <<Di  ,G0>r> 

<®2  >G<$>i)  {®2.G<5>2)  (<&2  >G<S>r) 

UO^GO,)  (<t>r,  G<t>2)  {<$>r.G<S>r)\ 


,G 


(Mb  +  vb^b) 
Mm  ”1”  Vm^rn 


(22) 


A 


diag 


(O^Oi)  1 
0 

0 


0 

(<F2>  ^2) 


0 


0 

0 

<<Dr,<I>r>' 


(23) 


The  polynomial  chaos  based  estimator  is  asymptotically  stable  if  the  eigenvalues  of  the  state 
transition  matrix  in  Equation  (20)  have  negative  real  parts.  For  the  base  excitation  model, 
stability  is  guaranteed  when  the  lower  bounds  on  the  parameter  prior  distributions  are  greater 
than  zero.  The  multivariate  polynomials  <t>k,k  =  1,2, ...  ,r  are  functions  of  the  Legendre 
polynomials  (Poularikas  (1999))  4>i(^m),i  =  0,1,  —,S  and  0;(£b)>7  —  0,1,  ...,S  as  follows: 

3 


k  —  i[  S  + 


+  7+1,  i  —  0,1, ... ,  S, 


J^D"  "  "U  '  2 

j  =  0,1,  ...,S  —  i 

The  inner  products  %B),  g(J;m,  ^b))  arc  defined  as  follows: 

1  1 


(24) 


(25) 


-1  -1 


Since  the  integrals  are  evaluated  over  the  event  space  of  the  random  variables,  the  inner 
product  (f  (^m,  ^g),  9  i.^m>  ^b))  ls  a  known  deterministic  quantity;  as  a  result,  Equation  (20)  is  a 
deterministic  set  of  state  equations  that  can  be  solved  numerically  to  determine  the  value  of  the 
state  vector  X.  Also,  since  the  multivariate  polynomials  /c  =  1,2,  ...,r  are  orthogonal  with 
respect  to  the  inner  product  in  Equation  (25),  the  matrix  Adiag  is  well  defined. 

Following  the  concepts  of  polynomial  chaos  theory,  the  states  xt  and  x2  of  Equation  (17)  are 
functions  of  the  random  variables  <f m  and  (B,  and  are  related  to  the  state  vector  X  of  Equation 
(20)  through  a  polynomial  chaos  expansion  approximation: 


■Xi 

x2. 


PX, 


P  := 


d>T 


0 


lxr 


0, 


o. 


9>r 


(26) 


'lxr  ^1 

This  approximation  becomes  exact  in  the  least  squares  sense  as  S  -»  00  (Ghanem  and  Spanos 
(1991)).  Substituting  Equations  (26)  and  (19)  into  (18)  results  in  the  polynomial  chaos 
approximation  output  y  which  is  a  function  of  the  unknown  variables  and  B . 

\-2{kf  +  kr)  -(jiB  +  vBfB)] 


y 


PX 


(27) 


L  flm  +  Vm^m  jlm  +  Vm^m  J 
The  output  y  is  not  an  output  trajectory  but  rather  an  output  process  or  family  of  trajectories. 
From  (27),  one  can  see  that  given  realizations  of  <f m  and  <fB,  y  collapses  to  an  output  trajectory. 
The  goal  of  recursive  maximum  likelihood  is  to  determine  the  realizations  of  and  ljB  that 
make  y  most  like  the  measured  output  y  in  some  sense.  This  is  done  by  selecting  and  tjB  so 
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that  they  optimize  the  likelihood  function  (assuming  zero  mean  Gaussian  measurement  noise 
with  variance  ok2)\ 

rfc  ,  ,  .  (  X'(y(tk)-y(tk'ZrrvtB)) 

I y)  °C  exp  j  -  >  - - 

v  k= 0  k 

Pence,  Lathy,  and  Stein  (2010  and  2011)  developed  techniques  for  recursively  calculating 
estimates  of  the  unknown  parameters  that  recursively  seek  to  optimize  the  likelihood  function 
(28).  The  key  limitation  of  the  polynomial  chaos  based  method  is  the  fact  that  the  polynomial 
chaos  approximation  is  only  exact  in  the  limiting  sense  as  S  -»  oo  which  is  not  numerically 
feasible.  The  authors  have  found,  as  a  rule  of  thumb,  that  S  >  6  is  a  good  approximation  for  this 
mass  estimation  problem.  A  general  procedure  to  determine  an  acceptable  value  for  S  is  to  start 
with  a  small  value  for  S  and  then  increase  the  value  until  the  change  in  the  resulting  estimates  is 
negligible. 


4  Regressor  Model  Based  Estimation  Methods 

An  alternative  approach  to  using  polynomial  chaos  based  estimation  is  to  use  regressor-model 
based  algorithms  such  as  recursive  least  squares  (RLS)  (Ioannou  and  Sun  (1996))  or  recursive 
total  least  squares  (RTLS)  (e.g.,  Kubus,  Kroger,  and  Wahl  (2008)).  Under  certain  Gaussian 
measurement  noise  assumptions,  and  if  the  states  xx  and  x2  of  Equations  (17)  are  measured 
explicitly,  RLS  methods  can  potentially  produce  unbiased  estimates.  Unfortunately,  full-state 
measurements  require  measuring  suspension  displacements  and  velocities  at  each  corner  of  the 
vehicle,  as  well  as  the  sprung  mass  acceleration.  This  paper  follows  Fathy,  Kang,  and  Stein's 
(2008)  method  in  using  pre-filtering  as  a  precursor  to  mass  estimation.  Such  pre-filtering  has 
two  very  attractive  advantages,  namely,  (a)  it  allows  the  estimation  process  to  focus  on  those 
frequencies  where  inertial  dynamics  are  more  visible,  and  (b)  it  makes  it  possible  to  estimate 
sprung  mass  using  sprung  and  unsprung  mass  accelerations,  without  any  need  to  use  additional 
sensors  for  displacement  and  velocity. 

RLS  and  RTLS  algorithms  rely  on  the  regressor  model  shown  in  Equation  (29)  and  cannot  be 
applied  directly  to  state-space  systems  such  as  the  system  modeled  by  Equations  (17)— (18). 

Yr  —  4>r  9r  (29) 

The  term  yR  is  the  measured  regressor  output ,  the  term  <pR  is  the  regressor  vector  which  is 

also  known/measured,  and  the  term  0R  is  the  unknown  parameter  vector.  Lor  vehicle  mass 
estimation,  the  regressor  model  can  be  derived  from  the  output  equation,  Equation  (18).  Then, 
the  regressor  output  yR  is  —  2{kf  +  /cr)x1,  the  regressor  vector  <pR  is  [7  x'z\T ,  and  the  unknown 
parameter  vector  6R  is  [m  BZ]T ,  i.e., 

—2 (kf  +  kr)x1  =  [y  x2]  [™]  (30) 

The  regressor  output  yR  and  regressor  vector  (pR  must  only  contain  known/measured 
variables,  but  the  states  xx  and  x2  are  not  measured.  To  address  this  problem,  this  paper  applies 
Laplace-domain  filtering  to  obtain  a  regressor  model  in  which  both  yR  and  <pR  only  contain 
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known/measured  values.  These  Laplace  domain  filters  act  as  pseudo-integrators  to  obtain 
estimates  of  velocity  and  displacement  from  measurements  of  suspension  acceleration. 
Assuming  zero  initial  conditions  for  the  states  x1  and  x2,  the  system  of  equations  given  in  (17)- 
(18)  can  be  represented  in  the  Laplace-domain  by  the  following  transfer  function: 

T(5)  _  Bzs  +  2(kf  +  kr) 

U  (s)  ms2  +  Bzs  +  2(kf  +  /cr) 

In  (31),  s  is  the  Laplace  operator,  Y (s)  is  the  Laplace  transform  of  the  sprung  mass 
accelerations  y,  and  U(s )  is  the  Laplace  transform  of  the  input  it.  Alternatively,  (31)  can  be 
written  as  follows: 

—2  (kf  +  /cr)(T(s)  —  U(s ))  =  ms2Y(s )  +  Bzs(Y(s )  —  U(s ))  (32) 

Dividing  both  sides  of  (32)  by  A(s),  results  in  the  following  filtered  regressor  model: 


— 2  (kf  +  kr) 


(TOO  -  t/(s)) 
AO) 


The  user-selected  denominator  A(s)  shou 


s2T(s)  s(Y(s)  —  U(s)) 


m- 

Bz. 


m  a o)  ■  (33) 

Id  be  a  polynomial  function  of  s  that  is  at  least 
second  order  to  ensure  that  the  transfer  functions  in  (33)  are  proper.  The  polynomial  A(s)  must 
also  have  complex  roots  with  negative  real  parts  to  guarantee  stability.  Reverting  back  to  the 
time  domain,  the  regressor  model  shown  in  Equation  (33)  becomes: 


-2(fc/  +  fcr)  Cy  -  u)  = 


L  ( 


Ia Tp) 


0-u)  [sj 


(34) 


Wmr 

The  term  (1/A(p)}(y  —  it)  represents  the  time  domain  signal  y  —  it  filtered  by  the  Laplace 
domain  transfer  function  (1/A(s)).  The  use  of  the  operator  p  instead  of  s  is  to  distinguish 
between  the  time  and  Laplace  domains. 

Equation  (34)  contains  only  known/measured  variables  in  the  regressor  output  yR  = 
— 2{kf  +  /cr){l/A(p)}(y  —  it)  and  regressor  vector  cj)R  =  [{p2/A(p)}y  (p/A(p)}(y  —  u)]T . 

The  only  unknown  variables  in  (34)  are  the  elements  of  the  parameter  vector  0R  —  [m  BZ]T . 
Therefore,  Equation  (34)  is  a  valid  model  for  these  regressor  model  based  algorithms.  However, 
two  significant  problems  arise  when  using  Equation  (34).  First,  both  the  regressor  output  yR  and 
regressor  vector  (pR  contain  measurements  of  both  y  and  it.  Thus,  the  noise  in  the  output  yR  and 
regressor  vector  <pR  is  correlated.  This  leads  to  biased  estimates  when  using  recursive  least 
squares  algorithms.  This  approach  was  explored  in  a  simulation  study  by  Pence,  Fathy,  and 
Stein  (2008a).  Total  least  squares  regression  can  potentially  lead  to  less  biased  estimates  (Moon 
and  Stirling  (2000)).  However,  a  second  problem  with  using  Equation  (34)  affects  both  least 
squares  and  total  least  squares  methods:  these  methods  are  both  sensitive  to  the  user-selected 
denominator  A  (.s’).  As  the  experimental  validation  section  will  show,  slightly  different  tunings  of 
A(s)  lead  to  significantly  different  estimates  of  the  unknown  parameters,  and  tunings  that  are 
optimal  for  one  terrain  are  not  necessarily  optimal  for  another. 
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5  Filtering  Methods 


Filtering  algorithms,  such  as  the  Extended  Kalman  Filter  (EKF),  Unscented  Kalman  Filter 
(UKF)  (Simon  (2006)),  particle  filter  (Ristic,  Masked  and  Gordon  (2004)),  etc.,  can  estimate 
both  the  states  and  unknown  parameters  of  state  space  systems.  Unlike  the  estimation  methods 
discussed  in  the  paragraphs  above  which  require  a  regressor  model,  the  filtering  algorithms  can 
be  applied  to  state  space  systems  such  as  the  base  excitation  system  of  Equations  (17)— (18). 
First,  however,  the  state  vector  must  be  augmented  to  include  estimates  of  the  unknown 
parameters.  This  augmented  set  of  state  and  output  equations  for  the  base  excitation  system  of 
Equations  (17)— (18)  is  given  as  follows: 


Xi 

*2 

*2 

2  (kf  +  kr)  x4 

- Xi - x2  —  it 

*3 

*3  *3 

X4_ 

O  O 

(35) 


2  (kf  +  kr)  x4 

y  - - xx - x2 

*3  *3 


(36) 


The  state  x3  is  an  estimate  of  the  unknown  value  of  the  vehicle  sprung  mass  m,  and  the  state 
x4  estimates  the  unknown  damping  term  Bz.  The  filtering  methods  discussed  above  can  be 
applied  to  this  augmented  system  of  Equations  (35)  and  (36).  This  augmented  system  has  a 
higher  dimension  and  is  significantly  more  nonlinear  than  the  (linear)  system  of  Equations  (17)- 
(18).  Similar  to  the  regressor  model  approaches,  the  filtering  methods  are  more  difficult  to  tune 
than  the  proposed  polynomial  chaos  approach.  The  filtering  methods  have  a  large  number  of 
user-defined  tuning  variables:  the  4x4  process  noise  covariance  matrix,  the  measurement  noise 
variance,  the  initial  conditions  for  the  4x4  estimate  covariance  matrix,  and  the  initial  conditions 
for  the  state  vector  (a  total  of  at  least  25  scalar  tuning  parameters).  Despite  being  difficult  to 
tune,  the  EKF  is  especially  attractive  because  of  its  low  computational  demand  compared  with 
the  polynomial  chaos  approach  and  even  the  UKF  and  particle  filter. 


6  Pitch  dynamics  investigation 

In  Section  2,  the  2nd  order  base-excitation  model  Equation  (17)  was  developed  assuming  the 
pitch  dynamics  had  negligible  contribution  to  the  estimation  of  the  sprung  mass.  In  an  effort  to 
investigate  the  validity  of  this  assumption,  a  maximum  likelihood  polynomial  chaos  mass 
estimator  was  developed  based  on  the  4th  order  model  in  (16)  that  includes  the  pitch  dynamics. 
The  performance  of  the  two  estimators  was  compared  via  a  numerical  simulation  of  a  nonlinear 
seven  degree-of-freedom  full  vehicle  model,  developed  by  Li,  Sandu,  and  Sandu  (2005),  that  was 
parameterized  similar  to  the  2001  Nissan  Altima,  GXE  sedan  (as  used  in  Section  7).  The 
simulation  induced  a  relatively  large  pitch  (about  8  degrees)  in  the  vehicle  by  driving  it  over  a  30 
centimeter  bump.  Figure  4  shows  the  relative  performance  of  the  two  estimators.  The  4th  order 
estimator  is  shown  in  (a)  and  converged  to  an  answer  within  less  than  1%  of  the  actual  vehicle’s 
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mass.  However,  the  performance  of  the  2nd  order  estimator,  shown  in  (b),  resulted  in  an  estimate 
with  about  a  12%  error. 

This  study  implies  that  inclusion  of  the  pitch  dynamics  in  the  estimator  can  increase  the 
accuracy  of  the  estimates  when  larger  pitch  is  experienced.  The  computation  cost  of  the  added 
pitch  dynamics  was  roughly  4x  that  of  the  simpler  2nd  order  estimator;  this  corresponds  with  the 
0(n  )  different  between  the  size  of  the  two  models.  Therefore,  an  additional  degree-of- freedom 
is  available  in  choosing  which  estimator  to  implement.  If  sufficient  computational  resources  are 
available  and  reasonably  large  pitches  are  expected  then  the  4th  order  estimator  is  recommended. 

Estimated  Sprung  Mass  (Maximum  Likelihood)  Estimated  Sprung  Mass  (Maximum  Likelihood) 


_ (a) _ _ (b) _ 

Figure  4:  Estimator  results  for  a  simulated  vehicle  driving  over  a  30cm  bump,  (a)  Shows  the  results  for  the  4th  order 
Maximum  Likelihood  estimator  that  includes  the  pitch  dynamics,  (b)  Shows  the  results  for  the  proposed  2nd  order 
Maximum  Likelihood  estimator  that  does  not  include  the  pitch  dynamics. 

The  remainder  of  this  paper  is  dedicated  to  the  experimental  validation  of  the  methods  of  this 
paper. 

7  Experimental  Setup 

The  vehicle  used  for  the  experiments  of  this  paper  was  a  2001  Nissan  Altima,  GXE  sedan. 
The  accelerometers  used  to  measure  both  the  sprung  and  unsprung  mass  accelerations  were 
Silicon  Designs,  Inc.  single-axis  accelerometers,  model  2210-005.  Unsprung  mass 
accelerometers  were  attached  to  the  suspension  struts  (Figure  5)  using  J.B.  Weld®  epoxy.  The 
sprung  mass  accelerometer  was  attached  to  the  base  of  the  cup  holder  between  the  front  seats 
(Figure  6)  using  Duro®  super-glue.  The  data  acquisitioning  system  was  a  National 
Instruments®  model  NI  USB-6221.  Data  samples  were  taken  at  a  500  Hz  or  1  kHz  rate,  filtered 
using  a  low  pass  filter  as  discussed  later,  and  then  downsampled  to  a  100  Hz  rate.  The 
accelerometers  were  connected  to  the  acquisitioning  system  in  differential  mode. 

A  Proform®  67650  Vehicle  Scale  System  was  used  to  measure  the  total  vehicle  weight 
(including  the  driver  and  acquisitioning  equipment  weight),  which  was  found  to  be  1460  kg 
(3220  lbs).  This  Vehicle  Scale  System  also  determined  the  Center  of  Gravity  (eg)  location  to  be 
a  distance  of  62%  of  the  wheelbase  length  (2.59  m)  forward  from  the  rear  axle.  The  tires  and 
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wheels  were  also  weighed  and  found  to  have  a  total  mass  of  70  kg  (150  lbs).  The  remaining 
mass  of  the  unsprung  components  (hub,  suspension,  axle,  etc.)  was  estimated  to  be  94  kg  (210 
lbs)  based  on  shipping  weight  calculations  published  on  an  internet  website:  Amazon.com™. 
These  unsprung  components  with  their  associated  shipping  weight  estimates  are  listed  in  Table  1. 
Based  on  the  measurements  of  the  tire  and  wheel  weight  and  the  estimated  value  of  the 
remaining  unsprung  mass,  the  total  sprung  mass  was  1296  kg  (2860  lbs). 


Figure  5:  Placement  of  the  unsprung  accelerometer. 


Figure  6:  Location  of  sprung  mass  accelerometer. 
Table  1:  Mass  of  unsprung  suspension  components. 


Unsprung  Suspension  Components  and  Corresponding  Shipping  Weights  (kg) 

Front  suspension  (one  corner) 

Rear  suspension  (one  corner) 

Hub  Assembly 

2.3  kg 

Hub  Assembly 

2.4  kg 

Vi  of  Front  Axle 

4.2  kg 

Knuckle 

6.4  kg 

Steering  Knuckle 

3.7  kg 

Strut  Assembly 

4.8  kg 

Strut  Assembly 

4.4  kg 

Brake  Drum 

4.4  kg 

Disc  Brake  Rotor 

6.3  kg 

Drum  Shoes 

1.0  kg 

Brake  Caliper  and  Pads 

2.9  kg 

Brake  Hardware 

0.8  kg 

Vi  of  Control  Arm 

2.2  kg 

Wheel  Cylinder 

0.5  kg 

Vi  of  Bottom  Control  Arm 

0.6  kg 

Total 

26  kg 

Total 

20.9  kg 

UNCLASSIFIED 


13 


UNCLASSIFIED 


The  stiffness  values  of  the  front  and  rear  suspensions  were  determined  respectively  to  be 
41,800  N/m  and  31,700  N/m  per  corner,  e.g.,  the  front-right  suspension  stiffness  was  41,800  N/m 
which  was  also  the  value  of  the  front-left  stiffness;  the  total  front  stiffness  was  83,600  N/m  and 
the  total  rear  stiffness  was  63,400  N/m.  These  values  were  determined  by  placing  a  known  mass 
of  53  kg  on  the  sprung  body  at  the  front  middle  and  then  rear  middle  of  the  vehicle;  at  each 
placement,  the  distance  between  a  fixed  point  on  each  tire  hub  and  a  fixed  point  directly  above 
on  the  body  of  the  vehicle  was  measured  using  a  ruler  with  millimeter  precision.  Each 
measurement  was  compared  with  its  unloaded  value  to  determine  the  respective  displacement  in 
the  suspension.  The  known  added  mass  and  suspension  displacement  values  were  assessed  with 
a  vertical  static  force  balance  to  determine  the  stiffness  values  of  the  suspensions.  Each 
suspension  displacement  for  the  two  different  placements  is  shown  in  Table  2. 


Table  2:  Suspension  displacement  under  a  53  kg  load. 


Placement 

location 

Front  Feft 
Displacement 

Front  Right 
Displacement 

Rear  Feft 
Displacement 

Rear  Right 
Displacement 

Front  Middle 

8.5  mm 

8.5  mm 

-3.0  mm 

-3.0  mm 

Rear  Middle 

-4.0  mm 

-4.0  mm 

13.5  mm 

13.5  mm 

Experimental  results  from  three  terrain  types  are  presented  in  this  paper.  The  three  terrains 
are  listed  in  order  from  smoothest  to  roughest:  smooth  pavement,  dirt/gravel  terrain,  and  rough 
pavement  with  potholes  and  cracks.  Figure  7  shows  frequency  responses  of  the  unsprung  and 
sprung  mass  accelerations  that  resulted  from  these  three  terrain  types.  The  target  speeds  were  40 
km/h  (25  mph)  for  the  smoothly  paved  terrain,  25  km/h  (15  mph)  for  the  dirt/gravel  terrain,  and 
40  km/h  (25  mph)  for  the  roughly  paved  terrain. 


(a)  Smooth  Hz  (b)  Dirt/Gravel  Hz  (c)  Rough 


Figure  7:  Spectrum  of  acceleration  magnitudes  for  (a)  smooth  pavement,  (b)  dirt/gravel  terrain,  (c)  rough  pavement. 

Figure  8  shows  a  frequency  response  plot  of  the  output  (sprung  mass)  acceleration  relative  to 
the  input  (averaged  unsprung  mass)  acceleration  for  the  dirt/gravel  terrain.  This  plot  was 
obtained  by  taking  the  discrete-Fourier  transform  of  the  output  divided  by  the  discrete-Fourier 
transform  of  the  input.  A  Bode  magnitude  plot  (solid  line)  of  the  base  excitation  model, 
Equation  (17),  with  m  —  1295  and  Bz  —  13000  is  also  shown  on  this  graph  as  a  reference.  This 
input-output  frequency  plot  shows  a  resonance  peak  between  0.9  -  4  Hz.  It  also  shows  increased 
magnitudes  for  frequencies  above  10  Hz.  These  large  magnitudes  in  the  higher  frequency  range 
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are  the  effect  of  disturbances  not  captured  by  the  base  excitation  model.  To  reduce  the  effect  of 
these  disturbances  on  the  parameter  estimates,  this  paper  applied  a  second  order,  low  pass  filter 


with  a  natural  frequency  of  5  Hz  and  a  damping  ratio  of  0.7  to  filter  the  sprung  and  unsprung 


acceleration  signals. 

The  recursive  estimation  algorithm  based  on  polynomial  chaos  theory  was  applied  to  the  data 
sets  from  the  three  terrains  discussed  above.  The  polynomial  order  was  set  to  be  S  —  6,  the 
bounds  on  the  mass  estimate  were  800  kg  <  m  <  2300  kg,  and  the  bounds  on  the  damping 
coefficient  were  4000  <  Bz  <  28000  Ns/m.  Figure  9  shows  the  convergence  of  the  polynomial 
chaos  algorithm  for  vehicle  sprung  mass  for  the  different  terrain  types.  For  each  terrain  type,  the 
algorithm  converged  to  within  10%  of  the  true  sprung  mass  value.  Assuming  convergence 
within  10%  error  to  be  acceptable,  these  results  verify  that  the  proposed  method  can  be  used 
successfully  for  vehicle  sprung  mass  estimation  on  both  smooth  and  rough  terrain. 


10 


10 


-3 


10 


-2 


10 


-1 


10 


10 


Hz 

Figure  8:  A  spectrum  of  input  to  output  magnitude  for  the  dirt/gravel  terrain.  Both  axes  are  shown  in  a  logarithmic  scale. 
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Figure  9:  Convergence  of  the  proposed  estimator  for  (a)  smooth  pavement,  (b)  dirt/gravel  terrain,  (c)  rough  pavement. 
The  dashed  line  is  the  true  mass  value,  and  the  dotted  lines  are  the  ±10%  error  lines. 
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Vary  nm  Varyvm 


Figure  10:  Sensitivity  of  the  polynomial  chaos  algorithm  to  varying  the  value  of  (A)  fim  G  {701, 1300, 1500, 1900}  with 
vm  =  700,  and  (B)  vm  G  {300,  700, 1400}  with  fim  =  1500.  The  dashed  line  is  the  true  mass  value,  and  the  dotted  lines 
are  the  ±10%  error  lines. 

The  polynomial  chaos  algorithm  requires  a  maximum  of  6  user  defined  tuning  parameters:  the 
prior  mean  values  /im  and  /tB,  the  prior  maximum  variance  values  vm  and  vB,  the  value  of  the 
polynomial  chaos  order  S,  and  the  measurement  noise  variance  ak.  If  ak  is  a  constant,  its  value 
does  not  affect  the  parameter  estimates.  In  that  case,  the  polynomial  chaos  approach  is  not 
sensitive  at  all  to  the  user  define  value  of  ok.  Figure  10  shows  the  sensitivity  of  the  polynomial 
chaos  algorithm  to  various  prior  mean  /im  and  variance  vm  values. 


8  Experimental  Comparison  with  Regressor  Methods 

Using  the  same  three  data  sets  discussed  above,  this  section  compares  two  regressor  methods, 
Recursive  Least  Squares  (RLS)  and  Recursive  Total  Least  Squares  (RTLS),  against  the 
polynomial  chaos  method  discussed  in  the  previous  section.  The  experimental  study  found  that 
the  regressor  methods  were  sensitive  to  the  tuning  of  the  parameters  in  the  user-selected 
denominator  A(s).  In  this  experiment,  the  denominator  was  chosen  to  be  A(s)  =  s2  +  2 (cos  + 
co2.  The  value  of  (  was  fixed  at  (  =  0.0024  for  RLS  and  (  =  0.707  for  RTLS,  and  the  value  of 
co  was  varied  between  6  and  12  rad/s  (0.95  -  1.9  Hz).  Figure  11  demonstrates  this  sensitivity  for 
the  RLS  and  RTLS  approaches  using  the  dirt/gravel  terrain  data. 


Figure  11:  Sensitivity  of  varying  to  for  (a)  RLS  and  (b)  RTLS  approaches  for  the  dirt/gravel  terrain.  The  dashed  line  is 
the  true  mass  value,  and  the  dotted  lines  are  the  ±10%  error  lines. 
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Figure  12:  The  convergence  of  the  (A)  RLS  and  (B)  RTLS  estimators  for  (a)  smooth  pavement,  (b)  dirt/gravel  terrain,  and 
(c)  rough  pavement.  The  dashed  line  is  the  true  mass  value,  and  the  dotted  lines  are  the  ±10%  error  lines. 

The  next  experimental  study  demonstrates  that  the  optimal  (in  the  mean  square  sense)  tuning 
of  A(s)  for  one  terrain  type  is  not  necessarily  optimal,  or  even  appropriate,  for  another  terrain 
type.  The  optimal  values  of  a>  and  (  for  the  rough  pavement  data  were  found  using  the 
optimization  algorithm  FMINSEARCH  in  Matlab®  version  2008a.  For  the  RLS  algorithm,  the 
optimal  values  were  oo  —  8.44  and  (  —  0.0024.  For  the  RTLS  algorithm,  the  optimal  values 
were  a>  —  15.4  and  (  —  0.105.  Figure  12  shows  the  convergence  of  the  algorithms  for  the  three 
terrain  types  using  these  estimator  tuning  values.  Although  these  values  are  optimal  for  the 
rough  pavement,  comparing  Figure  12  with  Figure  11  reveals  that  they  are  not  optimal  for  the 
dirt/gravel  terrain,  and  they  are  even  outside  of  the  ±10%  error  lines. 


9  Experimental  Comparison  with  Filtering  Methods 

This  section  applies  two  filtering  methods,  the  Extended  Kalman  Filter  (EKF)  and  the 
Unscented  Kalman  Filter  (UKF)  to  the  three  data  sets  discussed  above.  Similar  to  the  regressor 
methods,  the  filtering  methods  can  be  difficult  to  tune,  especially  considering  the  fact  mentioned 
above  that  they  each  require  at  least  25  scalar  values  to  be  tuned.  As  an  example,  Figure  13 
shows  the  sensitivity  of  these  tuning  methods  to  the  value  of  the  measurement  error  variance. 
Despite  the  difficulty,  the  EKF  and  UKF  algorithms  can  be  tuned  to  converge  to  estimates  that 
are  nearly  identical  to  those  produced  by  the  polynomial  chaos  maximum  likelihood  approach. 
This  is  shown  in  Figure  14. 


Figure  13:  Sensitivity  of  the  (A)  EKF  and  (B)  UKF  estimators  for  different  values  of  noise  variance  =  7  ■  10",  n  = 
1,3, 5, 8. 
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Figure  14:  Convergence  of  the  EKF,  UKF,  and  polynomial  chaos  (PC)  methods  for  (A)  smooth  pavement  and  (B) 
dirt/gravel  terrain. 

Since  the  EKF  algorithm  is  more  computationally  efficient  than  the  polynomial  chaos  and 
UKF  algorithms,  when  computational  resources  are  limited,  the  polynomial  chaos  approach  can 
help  to  tune  the  EKF  algorithm  offline,  and  then  the  EKF  algorithm  can  be  implemented  online. 


10  Conclusions 

This  paper  has  introduced  recursive  techniques  for  sprung  mass  estimation  of  vehicles  on 
rough  terrain.  It  proposed  a  base  excitation  model  that  treats  vertical  unsprung  accelerations, 
instead  of  the  terrain  profile  or  a  suspension  actuator,  as  the  input  to  the  estimation  model.  The 
experimental  results  demonstrated  that  the  proposed  method  is  viable  for  vehicle  sprung  mass 
estimation  for  vehicles  driving  on  rough  (and  relatively  smooth)  terrain.  No  prior  knowledge  of 
the  terrain  profile  was  required,  and  no  active/semi-active  suspension  was  required.  The 
polynomial  chaos  approach  was  compared  with  regressor  model  approaches  and  filtering 
approaches.  The  regressor  methods  are  concluded  as  not  appropriate  for  this  sprung  mass 
estimation  approach  because  of  their  sensitivity  to  tuning  parameters.  The  filtering  methods 
were  also  more  difficult  to  tune  than  the  polynomial  chaos  approach,  but  they  could  be  tuned  to 
converge  to  nearly  the  same  estimates  produced  via  the  polynomial  chaos  approach.  Based  on 
these  results,  the  polynomial  chaos  approach  is  recommended  when  sufficient  computational 
resources  are  available,  otherwise  the  extended  Kalman  filter  is  recommended. 
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